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I. INTRODUCTION 

When a large quantity of energy is injected on a very small volume over an instantaneous 
time, matter and energy diverge from the deposited volume. This phenomenon generates 
a very high pressure over the surrounding gas producing an explosion. When the pressure 
difference between the region of explosion and the matter at rest is sufficiently large, a 
strong explosion is produced and the discontinuity gives rise to a divergent shock wave. 
An explosive shock wave has a perfect spherical shape if the volume where the explosion 
took place is sufficiently small and there are no anisotropics or barriers in the surrounding 
medium. It was first Sedov m who formulated and solved analytically this problem [see also 

n u 

|8[ . The problem was also treated by Stanyukovich in his PhD dissertation [see |9|, for a full 
description of the general solution] and in less detail by Taylor 12j. Both authors derived the 
corresponding self-similar equations and obtained numerical solutions. As shown by Sedov 
8|, the self-similarity index of the problem is obtained by requiring the energy inside the 
shock wave to be a fundamental constant parameter. Under these circumstances, the famous 
Sedov-Taylor similarity index a = 2/5 is obtained J5|,|l3| for a similarity variable ^ = r/t", 
where r is the radial distance measured from the point of explosion and t represents time. 
The inverse physical situation is represented by a strong imploding spherical shock wave 



propagating into a uniform density medium. This was first inve stig ated by Guderley | 



and 



fully solved by Landau & Stanyukovich [see for example Isl. lol. Il3l|. In this case, a detailed 
construction of the self-similar solution was found by demanding the solution to pass through 
a singular point admitted by the hydrodynamical equations. The similarity variable in this 
case is given by r^ = r/(— t)" with a similarity index a = 0.68837 for a monoatomic gas. 
The first attempt to build similarity solutions for relativistic hydrodynamics was intro- 



duced by the pioneering work of Eltgroth 21, l3l| with the identification of a self-similar 
variable J := v/c, where v is the velocity of the fiow and c the speed of light. Soon 
after that, Blandford and McKee [l| solved the relativistic similarity problem for a very 
strong explosive shock wave. On their analysis, they introduced a similarity variable 
X = (1 — r/ct) [1 + 2(m + l)r^] , where F represents the Lorentz factor of the explosive 
shock front measured in the rest frame of the un-shocked gas and m is the similarity index. 
Their analysis represents a relativistic generalisation of the work carried out by Sedov- 
Taylor-Stanyukovich. 



In this article we construct self-similar solutions for the case of a very strong imploding 
shock wave converging to a centre. This constitutes a relativistic generalisation to the 
Guderley-Landau-Stanyukovich problem. In order to obtain a unique description of the 
motion, as it happens for the non-relativistic case, we force the solution to pass through an 
admissible singularity point, which in turn fixes the value of the similarity index of the flow. 

The present article is organised as follows. The next section describes briefly the main 
equations for relativistic hydrodynamics on a flat space-time and mentions the main results 
for shock waves on any inertial system of reference. Later in the article, we write down 
the self-similar equations behind the flow of a strong relativistic imploding shock wave. We 
calculate the similarity variable, the self-similar index of the problem and solve numerically 
the equations of motion. Finally, we describe the main properties of the solution. 

II. SHOCK WAVES IN RELATIVISTIC HYDRODYNAMICS 

The equations that govern the motion of a relativistic flow on a flat space-time are well 
described by Landau and Lifshitz [5|. First, the conservation of energy and momentum are 
represented by the condition that the divergence of the energy-momentum tensor T^'^ must 
vanish, i.e. 
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with 

T^""" = wu^'u'' - p7]^'\ (2) 

In here and in what follows Greek indices have values 0, 1, 2, 3 and Latin ones have values 
1, 2, 3. As usual, we use Einstein's summation convention over repeated indices. The 
coordinates x° are such that x° = (ct, r), where r is the spatial radius vector, t represents 
the time and c is the velocity of light. For a flat space-time, the metric tensor 7]^,^ is given 
by ?7oo = I7 ''711 = 1122 = "/^sa = — 1 and rjap = for a 7^ /?. The four-velocity u^ is 
defined by the relation m" = da;°/ds, where the interval ds is given by ds^ = rj^ydx^dx'^ . 
All thermodynamical quantities, in particular the pressure p, the enthalpy per unit volume 
w = e + p^ the internal energy density e and the temperature T are measured in the proper 



frame of the fluid. The internal energy density e constructed in this way contains the rest 
mass energy density pc^ and the "purely" thermodynamical energy density £, that is 

e = pc^ + e, (3) 

where p represents the proper mass density related to the proper particle number density n 
and average mass m per particle by the relation 

p = nm. (4) 

In the absence of sources and sinks of particles, the particle number density n satisfies 
the continuity equation 



5x° 



0, 



(5) 



where n" is the particle number density four-vector defined by n" := nw". 

For one dimensional flow, equation (^ can be rewritten as two equations, the equation 
of conservation of energy {dT^°'/dx°' = 0) and the equation of conservation of momentum 
(9T^"/9x" = 0) respectively as 
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(7) 



l-p\ 

Planar, cylindrical and spherical symmetric flows are described by these equations when 
k = 0, 1 and 2 respectively. 
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The relativistic theory of shock waves was first investigated by Taub 
well described by Landau and Lifshitz '5]. Traditionally the analysis is made on a system 
of reference in which the shock wave is at rest. In a more general approach, McKee and 
Colgate p wrote the necessary equations that describe a shock wave on any inertial system 
of reference. For the case of an imploding shock wave, it is convenient to choose a system of 
reference in which the un-shocked gas is at rest. In this particular system of reference the 
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jump conditions across a strong shock wave are simple and take the form [1|, |6|| 
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(8) 



- = 7 7 + 7, 9 

pi K — 1 K — 1 

r = ^ =■ (10) 

In the previous equations and in what follows, quantities without subindexes and those with 
subindex 1 label post-shocked and pre-shocked hydro dynamical quantities respectively. The 
Lorentz factor of the post-shocked flow and that of the shock wave are represented by 7 and 
r respectively. The corresponding values of the velocity measured in units of the speed of 
light are (3 and jSg. On both sides of the discontinuity, a Bondi- Wheeler equation of state 
with constant index k, given by 

p={K-l)e, (11) 

has been assumed. 

III. SELF SIMILAR EQUATIONS AND BOUNDARY CONDITIONS 

Let us consider a strong spherical shock wave converging to a centre. The pre-shocked 
material is assumed to have a constant density pi. We analyse the problem in a state such 
that the radius R of the shock wave is much smaller than the initial radius Rq it had when 
it was produced. Under these circumstances, as it happens for the non-relativistic case 
fl a fl, *Ke «ow is .a,e.y .dependent of tKe specie ...a. co.d.Uon. One ca. thus 
imagine a spherical piston that pushes the flow to the centre from very large distances in 
such a way that it produces a strong shock wave. That is, 

^ » ^. (12) 

n rii 

This expression combined with equations (jHI) and (P) for a gas with the equation of state (|TT|l 
implies that 

7 - 1 > 7 ^— 1 + 7 ^— , (13) 

K — I C^pi \ Ki — I C^PiJ 



where ki and k are the constant indexes for the pre-shocked and post-shocked gas respec- 
tively. The Lorentz factor 7 for the shocked gas that satisfies the previous inequahty is so 
large that the post-shocked flow velocity is comparable to that of light. 

In order to simplify the problem, let us assume that the particles that constitute the 
shocked-gas are themselves relativistic, so that they obey an ultrarelativistic equation of 
state with k = 4/3 in equation ()11|). That is, 

P = 36. (14) 

The shocked gas with this equation of state satisfies the following jump velocity condition 

7^ = ^^^ (15) 

according to equation (fTUI) . From this, it is clear that inequality ()13|) imposes the condition 
r^ ^ 1 for a strong shock wave. This justifies the restriction of our further analysis to 
zeroth order quantities in F^ and 7^. 

The density jump condition given by equation Q takes the form 

h = 47^ni, 

in which n represents the particle number density of the shocked gas as measured with 
respect to the pre-shocked gas, i.e. n = 'jn. In terms of the Lorentz factor F of the shock 
wave, this last equation becomes 

n = 2F^ni. (16) 

Finally, the energy jump condition (jH)) is given by 

n 2 

e = — 'juji = 2F uji. 
rii 

In terms of the pressure p, the previous equation is 

P=^rV, (17) 

where lji = pi(? for a cold pre-shocked gas and uji = Api for an ultrarelativistic one. The 
three conditions given by equations (fT3j). (fTBj) and (jTTj) represent post-shocked boundary 



conditions of the problem. These conditions are exactly the same ones obtained by Blandford 
and McKee [l| for the case of a relativistic explosive shock wave. 

Let us now write the equations of motion for the relativistic flow behind the imploding 
spherical shock wave. The symmetry of the problem demands that the energy conservation 
equation (jU)) and the momentum equation ((7j) take the form 



As a third equation we can either use the continuity equation, or the entropy conservation 
equation which is obtained from equations (fTHj). (jT^ and ^. We select here the last one, 
which for an ultrarelativistic gas is simply given by p| 

^ ^ =0. (20) 



dt Vn4/3 

Equations (fTHj) . (fT^ and (PUJ) together with the boundary conditions (fT3j). (fTBj) and (fT7|) 
completely determine the mathematical problem. 

IV. SELF-SIMILAR SOLUTION 

To find a self-similar variable for the problem, we recall that for the case of a relativistic 
strong explosion [l| , if i? represents the radius of the shock measured from the origin, then 
the ratio R/T^ is a characteristic parameter of the problem. With this ratio, Blandford and 
McKee [jj built a similarity variable for the problem. In the case of an imploding shock wave, 
none of the arguments in favour of the characteristic ratio R/T'^ are valid and instead we 
must construct the similarity variable by other m eans. The Appendix sh ows an alternative 



Blandford and McKee can be calculated. 



way in which the similarity variable obtained by 

The advantage of this approach is that we can generalise the method for the case we are 

dealing with in this article. 

The similarity variable for the case of a strong relativistic implosion is found as follows. 
The radius of the shock wave R converges to the origin as time increases, so that its velocity 



dR/dt must be negative. Thus, 

IdR 
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(21) 



In order to integrate this equation, we must know the time dependence of the Lorentz factor 
r. For the case of a strong relativistic explosion, this is given by equation ()A3j) . We can 
use a similar time power law dependence by taking into consideration the fact that the 
shock wave accelerates as it converges to a point. In other words, the Lorentz factor F must 
increase as time advances, that is 



A(-t)-™, 



(22) 



for times t < 0. The value of the constant factor A is determined by the specific initial 
conditions of the problem and m is the similarity index. The time interval t < has been 
chosen in complete analogy to the Newtonian case [cf . y, |l3[ . Under these circumstances, 
the time t = corresponds to a time in which the shock wave collapses to the origin, i.e. 
R = 0. 

Substitution of equation (j221) on relation (PT|) yields the integration 



R = c\t\ 



1 
1 



2(m + l)F2j- ^^^^ 

In exactly the same form in which we built the similarity variable for the case of a strong 
relativistic explosion (cf. Appendix), let us choose here the similarity variable r] as 



r/ =1 + 2(m + 1)F2 (1 - r/i?) 



:l + 2(m+l)F^ 



1 + 



1 



ct(l-l/2(m + l)F2)_ 



Thus, to order 0(F ^), this similarity variable takes the form 



V 



(l + ^) [l + 2(m+l)F2] 



(24) 



The problem we are concerned is such that the relevant region r > R corresponds to the 
values 1 > 1] > — oo of the similarity variable rj. 



For this kind of self-similar description, Buckingham's theorem of dimensional analysis 
motivated by the boundary conditions given by equations (fT3j) . (fTBj) and (fTTj) . demands 
the existence of three dimensionless hydrodynamical functions /(^), g{ri) and h{ri) given by 

p=luJiT'f{v), (25) 

n =2niT^h{ri), (26) 

f =^r^^(r^). (27) 

These dimensionless functions should be analytical and satisfy the following boundary con- 
dition at the shock wave 

f{r] = 1) = g{7^ = 1) = h{r] = 1) = 1. (28) 

The introduction of the similarity variable rj means that the hydrodynamical functions 
change their independent variables to F^ and rj, instead of r and t. The derivatives of r 
and t in terms of these new variables are thus given by 



^ m^J^+[l + {m + l){2T^-ri)]^, (29) 



dint SlnF^ '- ^ dr] 

ct|: = [l + 2(m + l)F^]|-. (30) 

Using these two equations, the total time derivative d/dt = d/dt + vd/dr takes the form 

m— — + (m + l) --r^ — . (31) 



dint 91nF^ \g J drj 

With the aid of equations (J25|l . (J26|l and (J27j) we can rewrite equations (fT8|l and (fT9|l in 
terms of F^ and rj. In fact, to 0(F^^) we obtain 



f q' 

(m + 1)(2 + gT])j^ + 4(m + 1)(2 - gr])^ = 3m, (32) 



^ [(m + 1) (8 - 2gri)] - 8(m + 1)^ 

/^ g 



4^ [(m + 1) (8 - 2^7/)] - 8(m + 1)^ = 2m - 8, (33) 



where the prime denotes derivative with respect to 77. These two equations can be rewritten 
in matrix form as 



-{2 + gr^) 2(2 -gr,) 
8 - 2gri -8 



r/fg 
g'/g' 



3m/{m + 1) 
(2m-8)/(m+l) 



which imphes, according to Cramer's rule, that 



(34) 



1 d/ 8(m - 1) + 5f?7(4 - m) 

fg dr] (m + 1) [—4 + Sgr] — g'^rf] 
1 dg 4 — 7m + grj{2 + m) 



(35) 
(36) 



g"^ dr] [m + 1) [—4 + 9>gri — g'^rf] 

Equations (jH^j) and (jHBjl are accompanied by the boundary conditions (f^ for firj) and girj). 

The unique value for the similarity index m and the uniqueness of the solution are de- 
termined by the following physical consideration, analogous to the one presented in the 
Newtonian case [cf . j^ . 

For a given time, the values of the energy density and velocity of a particular fluid 
element must decrease as the radius r increases. For the case of an ultrarelativistic gas, 
these quantities can be written in terms of the pressure p and the Lorentz factor 7 of the 
flow respectively. Thus, a necessary condition for the flow that guarantees the decrease of 
these quantities is given by 



dp 



<0, 



97' 



<0, 



^r > R. 



dr ' dr 

In terms of dimensionless quantities, the previous conditions are 



d/(r/) 



>0, 



dg{7]) 



>0, 



Vr/ < 1. 



(37) 



drj drj 

The integral functions /(?]) and g{ri) depend on the choice of the parameter m and so, 
a family of integral curves can be found for different self-similar flows. However, with the 
restrictions imposed by inequalities (j?7j) . it is possible to integrate the differential equations 
(jH^ and (jHBj) in a unique form. This happens because the conditions (jHTj) impose a single 
value for the similarity index m. 

The derivative dg{ri)/dri in equation ()36|) has a positive value at the surface of the shock 
wave, i.e. when rj = g = 1. In order to guarantee a positive value for all r/ < 1, g{rj) 
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must avoid extrema. In the same form, the continuity of the function g{ri) means that the 
denominator on the right hand side of equation (jHUj) should not vanish. To overcome these 
problems, we restrict the integral curve to pass through the singular point (77*, g{ri*)) for 
which the numerator and denominator on equation (j36p vanish simultaneously. 

The denominator on equation (|36p vanishes when the following condition is satisfied 



9V 



4 ± 2V3, (38) 



which represents a hyperbola on the grj plane. The numerator of the same equation has a 
null value on the hyperbola 

4 - 7m + g7]{2 + m) = 0, (39) 

The point {Ti*,g{ri*)) where both hyperbolas intersect is calculated by direct substitution 
of the value gr] on equation ()39|) . choosing the negative sign so that the number m remains 
positive. This fixes a unique value for the similarity index m given by 

m = I2V3 - 20 = 0.78460969. (40) 

If we now follow the same procedure for the derivative df{ri)/dri in equation (jH^j) . we 
obtain the same value for the similarity index m as the one given by the previous equation. 

The value of the similarity index given by equation ()40|1 guarantees that pressure and 
velocity gradients have negative values for the shocked gas. Using a 4th rank Runge-Kutta 
method we have calculated the integral curves for g{r]) and for /(//). These results are 
plotted in Fig. [U 

Let us now obtain a differential equation for the function h{r]). This is obtained by direct 
substitution of the differential equations (|HH|) and ^^ in equation (j^ . that is 

1 d/i _ 2 (8 - 9m) + 2gr] (5m - 6) - g^r]^ (m - 2) 
ghdrj [2 - gr]){m + 1) [-A + 8gr] - g^r]^] 

The integral curve /i(?7) of this equation shown in Figure Q was calculated with the already 

integrated values oi g{ri) and the unique value of m given by equation (j40j) . The number m is 

such that when the numerator on the right hand side of the previous equation vanishes, the 

numerator also does so, avoiding singularities and extrema in the range f] < 1. Furthermore, 

the derivative dh/drj evaluated on 77 = 1 has a negative value. Thus, the particle number 
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density n as presented by relation 
of time t, i.e. 



grows monotonically as r increases for fixed values 



dh{ri) 
dr] 



<0. 



(42) 



This peculiar situation also occurs for the non-relativistic imploding shock wave and will be 
discussed in more detail in the next section. 

With the integral functions /(r/), g{ri) and h{ri) calculated in this section it is then possible 
to calculate the dimensional hydrodynamical quantities using equations (j^ . (PB|) and (fTfj) . 

V. DISCUSSION 

The self-similar solution found in this article shows that some of the values of the hydro- 
dynamical quantities decrease in a given fluid element as it moves behind the shock wave. 
For example, the variation of the pressure in time as a fluid element moves behind the shock 
wave satisfies the following inequality 



dp 
dt 



1 dp 
|t|dlnt' 

1 r d 

t 



dlnr^ 



+ {m + l) {--ri) — }p, 



< - 



< - 



dr] 

l_dl 
gf drj 



—mp + p{m + 1) (2 — grj) 

— {m + l)p + p (m + 1) (2 — grj) 

p{m+l){l- gr]) ——- 

igf OV 



J_df 
gf dr] 



(43) 



The solution obtained in the previous section satisfies the condition grj < 1, for r/ < 1 
(cf. Figure Q). From this and the condition in equation ()37|). it follows that the total time 
derivative dp/dt has a negative value for any fluid element behind the shock wave. Using 
similar arguments it is easy to prove that, as the shocked gas particles move away from 
the shock wave, the Lorentz factor 7 decreases and the particle number density n increases. 
These results combined with the fact that n = jn imply that the proper particle number 
density n increases as the fluid particles move away from the shock wave, i.e. dn/dt > 0. 
This result is exactly what is obtained in the non-relativistic case [l3| . The particle number 
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density increases behind the shock and at sufficiently large distances behind the shock it 
must converge to a finite value. 

Se.f-s„„,.a. solut,™. of the second tvpe (see e., ZePdovioh and Rai.. Q) sucK a. 
the one found in this article satisfy two basic properties. The first one is the fact that the 
integral curve passes through a singular point (for a specific value of the similarity index). 
The second one is the existence of a curve r{t) on the r-t plane which corresponds to the 
singular point r/* and is itself a C_ characteristic bounding the region of influence. In order 
to verify this last statement for the self-similar solution found above, we proceed as follows. 
Let r]* = const, correspond to a value of the similarity variable evaluated at the singular 
point. Using equations (^^ and (^^ it is then straightforward to calculate the differential 
equation satisfied by the curve r^{t) evaluated at the singular point: 

at 0(r^^). On the other hand, at the singular point, equation ()36|) implies that the product 
r]^g^ = 2/ (2 + -\/3) for the case < r^ < 1 which corresponds to the light causal region 
R < r < c\t\, according to equations (j^H]) and (J21I). With this, equation (jl^ can be 
rewritten as 

dr* c ,_> 

-— = -c+ — ^r — . (45) 

dt 2(2 + v^)7*2 ^ ^ 

The curve r{t) that describes the trajectory of any C_ characteristic on the r-t plane is 
described by the following differential equation 



of any ( 

m 



dr a — 3 ,,„. 

dt 1 — a (3 

where a and /3 are the sound speed and the velocity of the fiow measured in units of the 

speed of light c respectively. For the self-similar solution we are dealing with in this article, 

/3 ~ — 1 + 1/27^. This implies that at 0(7^^), equation (|46|) takes the following form 

dr _ c (1 - a) 

dt ~~''+272(l + a)- ^'^^> 

Now, using the equation of state p4|l . it follows that the value of the velocity of sound in 
units of the speed of light is given by a = 1/v^. So, equation ^7\ has exactly the same form 
as equation (|^ . This means that the curve r^{t) that corresponds to the singular value 77* 
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is itself a C_ characteristic. In what follows we label this particular characteristic as C* . 
From the above results, an important condition about the causality of the flow behind the 
imploding shock wave is obtained. Since the curve r^,(t) that corresponds to the singular 
point r/^, does not intersect the imploding shock wave (except at the origin of the r-t plane) 
and because the C_ characteristics do not intersect one another, there is a region near the 
shock wave for which all C_ characteristics intersect the shock wave. Beyond the curve r^{t), 
the C_ characteristics do not intersect the imploding shock wave. This means that the flow 
farther away from the singular C*_ characteristic does not have any causal influence on the 
shock wave. These physical consequences on the causality are exactly the same as the ones 
that occur for the non-relativistic imploding shock wave. 

The energy contained inside a gas shell that moves in a self-similar manner can be 
calculated as follows. First, let us write down the value for the energy E contained on a 
shell of shocked gas. For an ultrarelativistic gas, the energy density T^^ measured on the 
system of reference in which the pre-shocked gas is at rest is given by p(47^ ~ !)• Thus, to 
0(7"^) inside the shell limited by the radius -Ro(^) and Ri{t) the energy is given by 

E{Ro{t), Ri{t),t) = lQn p-f^"^ dr. (48) 

Jro 

When the limits of integration are taken from the radii r = i? to r* = 
c\t\ (1 — ri^/[l + 2(m + l)r^]), i.e. the self-similar region of the flow, we obtain 

E =8..,r^ \tf (1 + 2(,n + l)r=)- f f, (1 - ^^,(,„\^)p )' *'• 

^MuJiT^\t\^-^{r],). (49) 

In these relations, the quantity 2(77=,,) is a constant and all quantities were approximated to 
0(r~^). Substitution of equation (^^ in equation (jiUj) yields 

E = 87rcuiA|t|^~"'S(r7,). (50) 

This equation, together with the value of the similarity index given by equation ()40p . means 
that the energy contained inside the shell is proportional to a positive power of time. In other 
words, the energy inside this self-similar shell tends to zero as the imploding shock wave 
collapses to the origin, which occurs when t ^ 0. As the shock wave converges to a point, 
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energy becomes concentrated right behind the shock front. However the dimensions of the 
self-similar region decrease and the previous analysis shows that the energy concentrated 
within this region also decreases. 

The constant A was introduced in equation (j22|) . Its value is determined by the initial 
conditions that generate the shock wave, in complete analogy with the Newtonian case (cf. 
|l3[). Thus, through a specified value of this constant, different flows behind the imploding 
shock wave are generated for different values of the initial pressure discontinuity. FiguresElto 
13] show examples of numerical integrations of the equations of motion. These solutions were 
calculated in natural units such that c = ks = 1, with fee representing Boltzmann's constant. 
The pre-shocked uniform medium was normalised in order to have a particle number density 
ni = 1, initial pressure pi = 1 and initial enthalpy per unit volume Wi = mnic^ = 1. 
These assumptions fix uniquely the units of mass, length and time for the specific problem. 
The initial pressure discontinuity that gave rise to the shock was assumed to have a value 
P2 = 100 at a distance Ro = 100. The approximations made in this article break down for 
values 7 < 10. As discussed before, the pressure and the Lorentz factor behind the flow 
decrease with increasing radius. The particle number density seems to grow without limit 
behind the shock wave. This unphysical behaviour occurs because the solutions discussed 
above are valid only for 0(7~^), and the region where the proper particle number density 
converge does not satisfy this approximation. Only a full self-similar solution of the problem, 
with no first order approximations in 7"^, will show that behaviour. 

For the particular case mentioned in this article, we have presented the last stages before 
the collapse of an imploding shock wave, which was generated by a very strong external 
pressure and where the complete treatment requires a relativistic analysis. At these final 
stages of the implosion, the flow pattern is not too sensitive of the specific initial conditions 
of the problem. 

Finally, we mention the regime under which the discussed solution remains valid. This 
is simply determined by the value of the post-shocked flow velocity. For sufficiently large 
distances behind the shock wave, the velocity of the gas f3 in units of the speed of light 
reduces so that its Lorentz factor 7 ~ 1. In this velocity regime, the approximations made 
on this article are no longer valid and the relativistic self-similar solution has no meaning. 



15 



Ackno^vledgments 

We gratefully thank the anonymous referee who pointed out the importance of the singu- 
lar C*_ characteristic related to the nature of the self-similar flow of the second type discussed 
in this article. J. C. Hidalgo acknowledges financial support from Fundacion UNAM and 
CONACyT (179026). S. Mendoza gratefully acknowledges financial support from CONA- 
CyT (41443) and DGAPA (INI 19203-3). 

APPENDIX A: RELATIVISTIC SIMILARITY VARIABLE FOR A STRONG EX- 
PLOSION 

Let us consider a relativistic strong explosion, diverging from the origin of coordinates. 
As mentioned in the article, similarity solutions were found by Blandford and McKee [l[. 
Let us show in here how to obtain their similarity variable with an alternative argument 
that can be extended naturally for the case of a strong relativistic implosion. 

The rate of change of the radius R{t) of the shock wave in units of the speed of light is 
given by 

IdR / 1\'" 

;dr = (^-pj ■ '*^' 

where F represents the Lorentz factor of the shock wave. A Taylor expansion of the right 
hand side of this equation to 0(F"^) gives 

Let us consider that F^ is a power law function of time only, so that 

r^ = ^. (A3) 

in which A is an unknown constant (but can be found from the specific initial conditions of 
the problem) and the constant m represents the similarity index. Substitution of equation 
i3|) on ()A2|) yields the integral 



« = ^|i-;;7TT2rjl- (*^' 
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The product ct has dimensions of length. Thus, we can write the dimensionless ratio 

R _ (b 



r ' 212 (m+1)' 

where we have introduced a new dimensionless quantity (p. To first order approximation, 
the inverse of the previous relation is given by 



r 



R ^ ' 2r2(m + l)' 
so that the dimensionless variable 6 is 



(l-^)2(m + ir, (A5) 



and vanishes when r = R. 

Self-similar motion means that all hydrodynamical functions can be written as the prod- 
uct of two dimensionless functions n(t) and 7r{r/R) J8|, |l3[, which is satisfied by the param- 
eter on equation ()A5|) . To avoid a null value of this variable at the shock radius, it is more 
convenient to take as similarity variable the function % := + 1, and not 0, i.e. 



X=l+^1--J2{m + 1)T\ (A6) 

The range of this similarity variable should be taken for x -^ 1- The value x = 1 corresponds 
to the shock wave. As described by Blandford and McKee fl!], the hydrodynamical equations 
simplify to great extent when this variable is used. Substitution of equations ()A3|) and ()A4|) 
on ()A6|) yields an expression for x at 0(r~2) given by 

X=(l-^)[l + 2(m+l)r2]. (A7) 
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FIG. 1: The figure shows, from bottom to top, the integral curves of the dimensionless pressure 
/(r/), squared Lorentz factor g{ri) and particle number density /i(r/) as functions of the similarity 
variable rj. 
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FIG. 2: The figure shows pressure p profiles behind a strong relativistic imploding shock wave. 
Natural units for the plot were used. The ambient medium was assummed to be of uniform density, 
with a particle number density ni = 1, a pressure pi = 1 and an enthalpy per unit volume wi = 1. 
The imploding shock wave was considered to be formed at a radius ro = 100, where the initial 
pressure discontinuity p2 = 100. From left to right, the figure shows the pressure distribution as 
a function of radius r for time values of 4, 5, 7 and 10 before the collapse of the shock wave. The 
envelope for the different values of the pressure just behind the shock wave is also shown. The 
pressure grows without limit as the shock wave converges to the origin. 
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FIG. 3: Particle number density n profiles for the conditions mentioned in Figure EJ are shown 
in the figure, as a function of the radial distance r. From left to right, the plotted curves show 
the particle number density distribution for times of 4, 5, 7 and 10 previous to the collapse of the 
shock wave respectively. The envelope curve shows the particle number density for the gas that 
flows just behind the shock wave. The particle number density n tends to infinity as the shock 
wave reaches the origin. Note that the particle number density profiles appear to grow without 
limit for sufficiently large radius behind the shock wave. This is a false result due to the fact that 
the Lorentz factor of the flow 7 is such that 7^ < 10 and the assumptions made on this article are 
no longer valid on this regime (see text for more information on this). 
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FIG. 4: Different Lorentz factor 7 profiles as function of the radial distance r are shown in the 
figure. The profiles are calculated for the shocked gas behind a strong imploding relativistic shock 
wave travelling through a constant gas medium with the same conditions as the one shown in 
Figure 121 This plot shows, from left to right, Lorentz factor values corresponding to times of 4, 5 
7 and 10 before the collapse of the imploding shock wave. The envelope curve shows the values of 
the Lorentz factor for particles that have just crossed the strong shock wave. As the plot shows, 
the Lorentz factor diverges as the shock wave collapses to the origin. Note that the solution breaks 
down for sufficiently small values of 7 such that the 0(7"^) approximation is no longer valid. 
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